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ABSTRACT 

A  quasi- three-dimensional  rotor  viscous  code  is  used  to 
predict  high  subsonic  flow  through  an  annular  cascade  of  turbine 
blades.  The  well  known  Baldwin- Lomax  turbulence  model  is  used  in 
the  program.  An  attempt  was  made  to  implement  a  new  turbulence 
model,  based  on  renormalization  group  theory  in  the  program.  This 
was  done  to  improve  the  prediction  of  the  boundary  layer  transition 
on  the  blade  surfaces  .  and  subsequent  wake  development.  The 
comparison  of  these  two  turbulence  models  with  experimental  data 
are  presented.  Pressure,  velocity  ratio,  flow  angle  distributions 
and  downstream  wake  predictions  were  studied  using  results  from 
RVCQ3D  (Rotor  Viscous  Code  Quasi -Three -Dimensional)  code.  The 
computed  results  showed  good  agreement  with  experiment  when 
comparing  the  blade  surface  local  static  pressure  to  inlet  total 
pressure  ratio  at  the  midspan  position  of  the  annular  turbine 
cascade.  The  computational  approach  used  to  implement  the 
turbulence  model  is  also  described. 
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I .   INTRODUCTION 

A.   BACKGROUND 

The  term  renormalization  belongs  to  quantum  physics.  It 
arises  when  the  discrete  formulation  for  particles  is  extended 
to  include  the  case  of  continuous  fields.  Initially  this  is 
done  by  formulating  the  discrete  problem  on  a  lattice.  In 
turbulence  theory,  when  the  continuum  limit  is  taken,  and  the 
Kolmogorov  spectrum  hypotheses  are  applied  in  the  wavenumber 
space,  the  infrared  divergence  of  turbulence  theory  arises. 
Renormalization  is  a  standard  method  of  removing  the  resulting 
singularities  and  this  procedure  is  described  in  McComb 
[Ref .  1] . 

In  an  earlier  paper  [Ref.  2],  Yakhot  and  Orszag 
have  proposed  the  application  of  renormalization  group  theory 
to  turbulence  modelling.  In  the  RNG  approach,  the  Reynolds 
stresses  are  accounted  for  by  repeatedly  recasting  the 
equation  of  motion  with  the  smallest  spatial  scales 
represented  by  effective  larger  scale  motions.  Following  each 
renormalization,  the  spectrum  is  narrowed  by  removal  of  a  band 
of  the  highest  remaining  frequencies.  When  the  process  is 
complete,  a  set  of  averaged  equations  are  obtained  along  with 
equations  for  an  effective  turbulent  viscosity. 


Lund  [Ref.  3]  used  an  algebraic  mixing- length 
model  to  calculate  provisional  values  of  eddy  viscosity  at 
each  point  normal  to  the  wall.  These  provisional  values  are 
then  processed  by  the  RNG  equation  and  the  actual  values  of 
eddy  viscosity  to  be  used  in  the  mean- flow  calculation  are 
determined.  The  RNG  equation  is  such  that  if  the  provisional 
value  of  the  eddy  viscosity  is  below  a  certain  threshold 
value,  which  is  determined  by  the  Heaviside  function,  the  eddy 
viscosity  will  be  zero.  As  the  provisional  values  of  eddy 
viscosity  increase  above  the  threshold,  the  RNG  equation 
returns  non-zero  values  of  eddy  viscosity  that  are  always  less 
than  the  provisional  values.  Since  the  provisional  values  of 
eddy  viscosity  increase  -with  Reynolds  number,  the  RNG  eddy 
viscosity  will  be  zero  for  low  Reynolds  numbers,  increase  for 
moderate  Reynolds  numbers,  and  finally  asymptote  to  a  fully 
turbulent  value  for  large  Reynolds  numbers.  This  type  of 
behavior  should  make  the  RNG  approach  well  suited  for 
transition  calculations.  The  heart  of  the  algebraic  RNG  model 
is  a  ramp  function  that  switches  on  when  the  provisional  value 
of  eddy  viscosity  is  above  the  threshold.  If  the  ramp  function 
is  on,  the  RNG  equation  takes  the  form  of  a  cubic  equation  to 
be  solved  for  the  eddy  viscosity.  Lund  also  mentioned  a 
transformation  of  the  cubic  equation  to  a  quartic  equation  and 
that  this  form  of  the  RNG  model  would  possibly  be  more  robust. 

Kirtley  [Ref.  4]  used  a  similar  quartic  equation 
to  study  complex  three-dimensional  turbomachinery  flows.  As 


the  RNG  equation  was  expressed  as  a  function  of  mean  strain 
rate,  the  model  was  then  referred  to  as  " energy- based " .  In  the 
procedure  of  solving  the  eddy  viscosity,  the  shear  stress  was 
introduced  into  the  quartic  equation  and  a  modified  mixing 
length  was  applied  in  the  wake  region.  Newton's  method  was 
used  to  find  the  roots  of  the  RNG  equation,  with  the  shear 
stress  determined  from  a  previous  value  of  eddy  viscosity. 

B .   PURPOSE 

The  aim  of  this  thesis  is  to  improve  the  viscous  flow 
analysis  in  turbomachinery  blade  rows.  The  Navier-Stokes 
equations  need  to  be  solved  and  these  are  simplified  by 
Reynolds  averaging  which  gives  rise  to  the  Reynolds  stress 
tensor  in  the  momentum  equations.  Turbulence  models  are  needed 
to  account  for  the  Reynolds  stress  tensor,  and  these  represent 
the  transport  of  momentum  due  to  turbulent  fluctuations.  Good 
turbulence  models  have  been  developed  for  relatively  narrow 
ranges  of  flows  with  excellent  accuracy,  but  they  can  not  be 
applied  with  confidence  to  complex  flows. 

For  computational  efficiency  and  code  simplicity,  the 
Baldwin- Lomax  algebraic  turbulence  model  works  very  well  for 
many  two-dimensional  flows.  This  model  calculates  the 
viscosity  induced  by  the  vorticity.  However,  it  contains  many 
empirically  determined  constants,  which  are  determined  from 
simple  two-dimensional  shear- f low  experiments.  These  constants 


are  the  near- wall  damping  of  turbulence,  transition  location, 
and  intermittency.  In  the  analysis  of  turbine  flow  fields,  the 
transition  region  of  the  boundary  layer  occupies  most  of  the 
suction- side  surface.  The  predictions  of  skin  friction  and 
heat  transfer  are  important  in  turbine  blades,  therefore,  less 
reliance  on  empiricism  is  desired. 

The  RNG  formalism,  cannot  as  yet  account  for  free- stream- 
turbulence  intensity  effects.  The  elimination  of 
inf initesimally  thin  high-wavenumber  bands  is  repeated  over 
the  renormalized  spectrum,  giving  rise  to  recurrence  relations 
which  link  the  low- Reynolds -number  regions  of  the  turbulence 
with  the  high.  Some  constants  appear  in  these  relations  but 
they  are  determined  analytically  from  assumptions  about  the 
shape  of  the  energy  spectrum.  The  eddy  viscosity  is  a  function 
of  the  mean  strain  rate. 

The  object  of  this  study  is  to  apply  the  RNG  turbulence 
model  into  the  Rotor  Viscous  Code  Quasi -3-D  (RVCQ3D)  [Ref .5 
and  Ref.  6]  and  the  comparison  of  the  results  to  those 
computed  with  the  Baldwin- Lomax  turbulence  model.  Quasi -three- 
dimensional  flow  is  computed  by  numerically  solving  the  thin- 
layer  Navier-Stokes  equations  using  a  multi-stage  Runge-Kutta 
scheme  with  a  spatially  varying  time  step  and  implicit 
residual  smoothing. 

This  study  entailed  the  computation  of  two  cascade  test 
cases.  The  first  test  case  is  an  annular  cascade  of  turbine 
stator  blades  which  was  experimentally  measured  by  Goldman  and 


Seasholtz  [Ref.7]  with  a  one -component  laser  Doppler 
velocimeter.  They  mapped  out  the  flow  at  various  radial  and 
axial  stations  by  performing  pitchwise  blade-to-blade  and  wake 
surveys.  RVCQ3D  was  used  to  compute  the  midspan  location  and 
comparisons  were  made  with  the  experimental  measurements. 

The  second  test  case  is  a  subsonic  cascade  of  turbine 
rotor  blades,  experimentally  investigated  by  Hoheisel 
[Ref.8].  The  wake  traverse  measurements  were  performed 
with  a  wedge  type  probe,  and  the  boundary  layer  measurements 
were  carried  out  on  the  blade  suction  surface  with  a  flattened 
Pitot  probe.  This  test  case  was  only  computed  with  the 
Baldwin- Lomax  turbulence  model  and  is  presented  as  an 
Appendix. 


II.   FORMULATION 

A.   GOVERNING  EQUATIONS 

The  mass -averaged,  compressible  form  of  the  Navier- Stokes 
equations  theoretically  model  turbulent  flow.  The  equations 
are  solved  numerically  on  a  body- conforming  grid  system,  so 
these  equations  are  transformed  from  the  Cartesian  system 
(x,y)  to  a  body- fitted  coordinate  system  (£,tj)  .  The  thin- layer 
approximation  which  neglects  the  viscous -diffusion  terms 
parallel  to  the  body  surface  is  applied  after  the 
transformation. 

Liu,  Sockol,  and  -  Prahl  [Ref .9]  defined  the 
nondimensionalized  reference  quantities  as  follows: 

C* 

—    U  —  V  —       p 

U  =  —  ,        V   =  — ,   p  = 


C       (p0e-J 

T=J1,   e=  *— ,  (1) 

To  (p0c<2) 

P  =  -fi-,  F-JS-. 

Po         Ho 


N  Y+l        Ho  * 


The  nondimensionalized,  conservative  form  of  the  two- 
dimensional,  time -dependent,  mass -averaged,  thin- layer  Navier- 
Stokes  equations,  from  Chima  [Ref.  5],  are  as  follows 
(omitting  the  overbar) : 

Ti iJ'^  *J~k ( J"lF')  *Ji{  (J"1F^  -J5 J-k  lJ'lF^  '-  F>    (2) 

where  J  is  the  transformation  Jacobian  given  by 


J  = 


dx 

dx 

dy 

dy 

an 

(3) 


Equation  (2)  can  be  rewritten  as  follows: 


.A  ip )  +_i.  (p )  +_i_  ip  \  -_LjL  (p )  =p 


(4) 


to  which   the   numerical   scheme   is  applied 
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F2    = 


pv 
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0 

0 
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e     1 
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where 


and 


mi   =    (-|lx+lj)  (H+Hc) 


m2     =     T"TlxT1y(^+^C) 


»4     =      W(^^) 


(6) 


e  =  p    [  Cvr  +  l(u2  +  v2)]  (7) 


p  =   (y-1)  [  e  -  ^p(u2  ♦  v2)]  (8) 


(9) 


B.   TURBULENCE  MODEL 

1.   BALDWIN -LOMAX  MODEL 

This  model  is  a  two- layer  algebraic  eddy  viscosity 
model  proposed  by  Baldwin  and  Lomax  [Ref.10],  and  has 
been  extensively  used  for  calculating  turbulent  flows.  The 
description  of  the  model  follows. 


The  eddy  viscosity  is  calculated  using 


J  'rt'  inner  J     J crossover 

^  '  I  (l*e)  outer        ^     Y>Vc 


(10) 


crossover 


where  yerooover  is  given  by  the  minimum  value  of  y  for  which    (/xt! 

=   '  Mt '  outer  • 

The  inner  and  outer  values  are  given  by 


-2- 

A*\  \  2 


d*t)iM.r  =  P(Ky(i-e  A  ))2  » 


(11) 


where 


'r^c^outar  "  ^-outerP^WAKE^KLEB 


(12) 


y 


(13) 


A+  -  26,  C^  =  0.02  69,  and  FWAKE  is  given  by 

*Ws  =  min^max^nax  ,0.25  Y^V^2 / F^) 

The  quantities  y,,^  and  F,^  are  found  from 


(14) 


F(y)    =  y|o)|(l-e  A*) 


(15) 


where  y^  is  the  y  value  for  which  F(y)  is  a  maximum. 


Fkleb   is   given  by 


0.3  y 


Fkleb  =    [1+5. SC^-^)6]"1  (16) 


The  term  vDIF  is  given  by 


Except  for  wakes,  the  second  term  is  zero. 

A  mixing  length  is  introduced  in  the  inner  viscosity 
equation,  as  follows; 

2=Ky(l-e"0  U8) 

and  the  w  term  is  obtained  from  the  vorticity,  as: 

o>  =  1  (vorticity)    =  A(|Z-|H)  (19) 

2  2     dx     dy 

Since  the  coordinates  are  transformed  from  (x,y)  to 
(£,*;),  the  vorticity  in  the  equation  above  becomes 

b-*-?i^S^€>-<-^^S>]      (20) 

From  this  vorticity  form,  the  eddy  viscosity  can  be 
calculated  by  combining  the  inner  and  outer  turbulent 
viscosities . 
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2.   RNG  MODEL 

In  Kirtley  [Ref.  4],  the  following  relation  for  the 
renormalized  viscosity  is  derived: 

v  =  v0[l*H(-±*   A/4-Cc)]3 
v  J 
0        1  (21) 

v0 

where  H(x)  is  the  Heaviside  function,  which  is  zero  for 
negative  x  and  equal  to  x  for  all  positive  x.  \  is  the 
wavenumber  corresponding  to  the  largest  length  scale  in  the 
inertial  range.  The  coefficient  a  ■  0.1186  and  was  derived 
[Ref.  2]  from  the  -5/3  power  law  for  the  turbulence  energy 
spectrum.  Cc  is  a  function  of  both  a  and  the  ultraviolet 
dissipation  range  cut-off  in  the  spectrum  of  turbulence;  the 
suggested  range  is  between  75  and  200.  The  dissipation  rate  of 
turbulence  is  determined  from  the  following  equations: 

e  =  v  s  (22) 


where 


5=   ( i  + 2-_£$. £) i 

dXj     dXt     3    1J  dxk     dXj 

ox  ay  ox     oy         3     ox     oy 


(23) 
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Using  the  same  transformation  techniques  as  in  the 
discussion  previously,  the  mean  strain  rate  squared  is 
transformed  from  (x,y)  to  (£,17)  as  follows: 


As  discussed  in  reference  2,  the  nonlinear 
renormalized  viscosity  equation  .can  be  written  in  a  cubic 
equation  form: 

v3-vSLf4+(Cc-l)v03=0  (25) 

Kirtley  [Ref.  4]  showed  the  locus  of  possible 
allowable  roots.  If  the  maximum  value  of  v  is  chosen  as  the 
proper  root,  a  jump  discontinuity  in  the  distribution  of  v  is 
yielded.  A  quartic  equation  transformed  from  equation  (25)  was 
used  by  Lund  [Ref.  3],  as  follows: 

v4  -  vv03  -  Hix^Lf  -   vCcv03)  =  0  (26) 

where 


1 
t'  =  vS 


/-««-2  (27) 
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and  a  preferable  length  scale  proposed  by  Kirtley  [Ref .  4]  is 
the  following: 


Lf=C^   tanh(-^)  (28) 


with  CM  ■  0.0845,  k  =  0.372,  and  n  is  the  normal  distance  from 
the  boundary,  and  5  is  the  thickness  of  the  boundary  layer.  An 
approximation  of  <5,  was  suggested  by  Stock  and  Haase 
[Ref. 11],  as  follows: 

5  =  1.936  ymax  (29) 

where  y,,,,,  is  obtained  from  the  Baldwin -Lomax  turbulence  model. 
The  length  scale  Lf  in  the  wake  region  of  turbomachinery  flows 
is  determined  from  a  correlation  derived  by  Raj  and 
Lakshminarayana  [Ref. 12]  as  follows: 

Le  =  mindcs  ,  CJd)  (30) 

where 


b   =  60+c  Cd2  1.35  (-+0.02)0-58  <31> 

0  c 


and  50'  is  the  average  of  the  pressure- side  and  suction- side 
boundary  layer  thicknesses  at  the  trailing  edge,  c  is  the 
chord,  and  Cw,  the  Gaussian  constant  for  the  wake  behind  a 
circular  cylinder,  is  equal  to  0.169.  The  equations  above  are 
the  only  empiricism  in  the  RNG  model. 

Thus,  the  viscosity  on  every  grid  point  in  the  flow 
field  can  be  calculated  from  the  solution  of  the  quartic 
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equation  as  follows.  First,  an  initial  guess  of  the  turbulent 
eddy  viscosity  is  set  equal  to  zero,  before  the  flow  solver 
starts  its  time  marching  iterative  procedure.  The  algorithm 
for  the  RNG  turbulence  model  is  then  as  follows: 

1.  Calculate  the  boundary  layer  thickness,  mixing 

length,  and  mean  strain  rate  on  every  grid  point. 

2.  Calculate  the  shear  stress  to  density  ratio  based 

on  the  effective  eddy  viscosity  at  the  n  time  step. 

3.  Evaluate  the  Heaviside  function. 

4.  Apply  Newton's  method  to  solve  the  quartic  equation 
for  the  effective  eddy  viscosity.  There  are  two 
posibilities: 

a.  If  the  argument  in  the  Heaviside  function  is 
positive,  use  that  value  for  the  Heaviside 
function. 

b.  If  the  argument  in  the  Heaviside  function  is 
negative,  then  set  the  Heaviside  function  to 
zero. 

5.  Update  the  eddy  viscosity  and  return  to  step  3, 

until  the  solution  finally  converges. 

6.  The  converged  solution  is  the  effective  eddy 
viscosity  at  the  n+1  time  step. 

7.  Finally  reevaluate  the  Heaviside  function  again  and 
a.  if  H>0,  the  n+1  time  step  effective  eddy 

viscosity  is  used. 
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b.  if  H<0,  the  n+l  time  step  effective  eddy 

viscosity  is  set  equal  to  the  laminar  viscosity. 
This  means  the  turbulent  eddy  viscosity  equals 
zero. 

8.  A  turbulent  eddy  viscosity- to-dynamic  viscosity 

ratio  is  then  calculated  for  RVCQ3D. 

9 .  Return  to  main  program. 

After  the  RNG  turbulence  subroutine  calculates  the 
turbulent  eddy  viscosity,  the  value  is  stored  in  an  array  and 
becomes  the  initial  guess  for  the  next  iteration  calculation. 
The  calculation  direction  is  from  the  leading  edge  toward  the 
trailing  edge  and  then  on  through  the  wake.  This  was  done  so 
that  the  rate  of  growth  of  the  boundary  layer  could  be 
monitored. 

The  RNG  model  mimics  the  transition  from  laminar  to 
turbulent  flow  within  a  boundary  layer.  As  the  Heaviside 
function  suppresses  deviations  from  laminar  viscosity  until 
the  length  scale  and  mean  strain  grows  to  values  above  the  Cc 
cutoff.  The  same  result  occurs  within  the  turbulent- to- laminar 
transition,  region.  Thus  the  RNG  model  should  be  able  to 
capture  the  viscous  sublayer  without  any  near -wall  damping, 
while  the  Baldwin- Lomax  model  requires  the  addition  of  a 
damping  term  to  describe  the  near- wall  behavior. 
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III.   SOLUTION  METHOD 

A.   GRID  GENERATION 

A  transformation  was  used  to  transform  a  non- rectangular 
grid  in  the  physical  plane  into  a  rectangular  uniformly- spaced 
grid  in  the  computational  plane. 

Two-dimensional  body- fitted  grids  were  generated  using  the 
GRAPE  code  (GRids  about  Airfoils  using  Poisson's  Equation). 
This  code  was  written  by  Sorenson,  primarily  for  isolated 
airfoils  [Ref.13].  A  modified  GRAPE  code 
[Ref.14]  was  used  for  the  generation  of  periodic  C- 
type  grids  for  turbomachinery  cascades.  For  viscous  flow,  the 
grids  are  typically  finer  normal  to  the  blade  than  in  the 
streamwise  direction. 

The  locations  of  the  intersection  of  grid  lines  with  the 
inner  and  outer  boundaries  were  controlled  to  give  minimum 
shear  in  these  regions.  It  is  desired  to  have  near- wall  normal 
grid  lines  at  the  blade  surfaces.  However,  because  of  the 
large  amount  of  turning,  normal  grid  lines  could  not  be 
generated  over  the  entire  blade  surface  and  still  maintain  the 
grid  line  periodicity  on  the  outer  boundary.  The  presence  of 
the  sheared  grid  lines  did  not  appear  to  adversely  effect  the 
flow  solutions. 
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B.   NUMERICAL  SCHEME 

The  code  employs  an  explicit,  finite-difference,  multi- 
stage, Runge-Kutta  algorithm.  The  governing  equations  are 
discretized  using  a  node-centered  finite-difference  scheme 
[Ref .  15  and  Ref .  16]  : 

qx   =  F{Qn)   -  cc.At  R(F{0n))  (32) 

q2   =  F{Qn)    -  a2At  R(qx)  (33) 

g3  =  F{0n)   -  a3At  R(q2)  (34) 

g4  =  F(Qn)    -  At  R(q3)  (35) 

Pf1*     "  <U  (36) 

where  a3  =  1/2,  a2  -  1/3,  a,  =  1/4,  and 

(37) 

Re      4i    4j- 

Second-  and  fourth-order  artificial -dissipation  terms  are 
added  to  stabilize  the  scheme.  A  spatially  varying  time  step 
and  implicit  residual  smoothing  are  used  to  accelerate 
convergence . 
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C.   BOUNDARY  CONDITIONS 

At  the  cascade  inlet,  absolute  total  temperature  T0, 
absolute  total  pressure  P0,  relative  velocity,  and  relative 
flow  angle  are  specified  as  constants.  For  subsonic  flows,  the 
inlet  conditions  are  updated  at  each  iteration  by 
extrapolating  the  upstream- running  Riemann  invariant  R"  based 
on  the  absolute  total  velocity  Q  to  the  inlet,  where 

R-   =  Q  -      2c,  (38) 

Y-l 

Then,  using  the  known  total  temperature  T0  and  the  isentropic 
relations,  the  total  velocity  magnitude  is  given  by 


(y-l)R-  *  I  4(y+l)  CpTQ   -  2  (y-l)  _(Jr)^       (39) 

(Y+D 

Individual  velocity  components  are  found  from  trigonometric 
relations  by  keeping  the  tangential  velocity  component  fixed, 
Pressure  and  density  are  found  from  isentropic  relations.  On 
the  blade  surfaces,  for  viscous  flow,  the  no- slip  velocity 
boundary  condition  is  used.  The  surface  pressures  are  computed 
using  the  normal  momentum  equation,  namely 

(r\xZx   +  T\yZ.y)pK   +  (r\x2   +  tiy2)pn  =  -p*7(Tixu5  +  r\yvK)  (40) 

where  U  =  V  =  0  on  the  surface  for  viscous  flow. 

Surface  densities  are  found  from  surface  pressures  and  a 
specified  wall  temperature.  Both  the  Baldwin- Lomax  and  RNG 
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turbulence  models  incorporated  Sutherland's  viscosity  law  to 
calculate  the  dynamic  viscosity. 

The  exit  static  pressure  is  determined  by  specifying  the 
exit  static-to-inlet  total  pressure  ratio.  All  other  flow 
quantities  are  extrapolated  at  the  exit  plane.  Periodicity  is 
enforced  on  the  "C"  grid  by  using  a  pseudo-point  along  the 
periodic  boundaries  on  which  corresponding  pitchwise  variables 
are  located. 
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IV.   RESULTS  AND  DISCUSSION 

Two  blade  profiles  that  were  part  of  the  full -annular  ring 
of  36  blades,  with  C-type  gridding  over  the  computational 
domain,  are  shown  in  Figure  1.  The  chord  length  of  the  actual 
vane  was  55.54  mm  long.  The  geometry  of  a  full  250  X  60  mesh 
size  (150  grid  lines  around  the  airfoil,  50  grid  lines  on  each 
side  of  the  center  line  in  the  wake  region)  is  shown  in  Figure 
Al  in  Appendix  A. 


■ 

1. 1 


Figure  1   The  geometry  of  Goldman's  annular  cascade 
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The  experimental  investigation  of  this  cascade  [Ref.  7] 
was  conducted  with  the  hub-static  to  inlet-total  pressure 
ratio  maintained  at  a  value  of  0.65.  This  corresponds  to  a 
mean- radius,  ideal,  exit,  critical  velocity  ratio  of  0.78.  In 
Chima  [Ref.  17],  a  boundary  condition  and  initial 
condition  input  file  for  the  RVCQ3D  code  was  provided.  The 
computational  results  in  the  next  section  were  calculated  with 
an  exit-static  pressure  to  inlet-total  pressure  ratio  of 
0.685.  The  data  of  streamtube  variation  (Appendix  A.)  were 
provided  by  Goldman  [Ref.  18]  . 
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A.   STREAMTUBE  VARIATION  AND  GRID  RESOLUTION 
1.   STREAMTUBE  VARIATION 

The  predicted  blade  surface  pressure  distribution, 
without  streamtube  variation,  did  not  agree  well  wi'th 
Goldman's  experimental  data  [Ref.  7].  In  order  to  obtain  a 
good  baseline  solution,  with  the  Baldwin- Lomax  turbulence 
model,  streamtube  variation  based  on  Goldman's  [Ref.  18] 
throughflow  calculations  was  introduced. 

The  effect  of  streamtube  variation  in  the  convergence 
rate  is  shown  in  Figure  2.  The  application  of  a  varying 
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streamtube  resulted  in  better  convergence  as  it  maintained  a 
near  constant  convergence  rate,  particularly  after  1000 
iterations,  and  finally  gave  an  order  of  magnitude  improvement 
on  the  convergence  after  3000  iterations. 

Streamtube  variation  had  a  significant  effect  on  the 
prediction  of  the  blade  surface  pressure  distribution  as  can 
be  seen  in  Figure  3 . 


Figure  3    Streamtube  variation  effects  on  pressure  ratio 
distribution 


Constant  streamtube  thickness  and  radial  location 
versus  axial  distance  through  the  cascade  resulted  in  a  poor 
prediction  of  static  pressure.  Not  only  were  the  levels  of 
static  pressure  qu-ite  different,  but  the  shape  of  the  suction 
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side  pressure  distribution  was  also  not  predicted  properly. 

Finally  the  wake  distribution  predicted  with  and 
without  streamtube  variation  is  shown  in  Figure  4.  Once  again 
a  significant  improvement  on  not  only  the  level  of  velocity  in 
the  wake  but  the  shape  of  the  wake  profile  is  obtained  with 
streamtube  variation.  The  same  streamtube  variation  study 
using  a  grid  size  of  97  X  31  is  shown  in  Appendix  B,  Figures 
Bl,  B2,  and  B3 . 
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2.  GRID  RESOLUTION 

In  order  to  determine  the  best  grid  resolution  for  the 
blade-to-blade  viscous  flow  analysis,  the  influence  of  the 
grid  size  on  the  flow  characteristics  was  studied.  The 
solutions  of  three  different  grids,  97  X  31,  200  X  31,  and  250 
X  60,  were  compared  for  the  analysis  of  Goldman's  annular 
turbine  cascade. 

All  of  the  three  different  sized  grids  were  run  to 
3000  iterations  with  the  same  boundary  conditions  and  a 
Courant  number  of  4.0.  The  convergence  history  is  shown  in 
Figure  5. 
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Initially  the  convergence  rate  of  the  coarse  grid  is 
better  than  the  two  finer  grids,  however  this  shows  signs  of 
instability  in  the  latter  part  of  the  convergence  history. 
Both  the  fine  grids  (200  X  31  and  250  X  60)  showed  similar 
convergence  histories  and  they  both  achieved  approximately 
three  orders  of  magnitude  convergence  after  3000  iterations. 

The  effect  of  grid  resolution  on  the  prediction  of  the 
midspan  blade  surface  pressure  distribution  is  shown  in  Figure 
6. 
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Grid  refinement  from  97  X  31  to  200  X  31  shows  an 
improvement  in  blade  surface  pressure  distribution  when 
compared  with  experiment,  particularly  on  the  suction  side  of 
the  blade.  Further  refinement  from  the  200  X  31  to  250  X  60 
showed  little  improvement  in  the  pressure  distribution,  hence 
the  solution  with  this  size  grid  was  felt  to  be  grid 
independent  and  subsequently  used  throughout  the 
investigation.  The  250  X  60  grid  solution  also  gave  the  best 
resolution  at  the  trailing  edge  as  the  "overshoot"  in  P./PO^ 
was  the  least  for  this  test  case. 

An  investigation  of  the  grid  resolution  was  also 
performed  on  the  wake  profile  prediction.  At  the  153.2  percent 
axial  chord  position,  a  comparison  of  experimental  data  with 
the  prediction  of  local  velocity  to  critical  velocity  ratio  is 
shown  in  Figure  7 . 

The  coarse  grid  clearly  produced  a  wake  profile  which 
was  highly  grid  dependent,  as  the  profile  was  not  smooth  at 
the  wake  center  line.  Grid  refinement  improved  this  as  well  as 
the  prediction  of  the  free-stream  velocity.  However,  the 
pressure  side  profile  shows  an  overshoot  which  is  not  evident 
in  the  experimental  data.  The  250  X  60  grid  produced  the  best 
shape  for  the  wake  profile,  but  the  deficit  in  the  wake  is 
overpredicted  and  the  wake  diffusion  or  growth  is 
underpredicted.  The  discrepancy  in  the  free- stream  velocity  is 
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attributed  to  the  inability  of  the  streamtube  contraction  to 
account  for  the  three-dimensional  effects  in  the  annular 
cascade. 
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B.   B-L  MODEL  AND  RNG  MODEL 

In  order  to  test  the  performance  of  the  RNG-based 
turbulence  model  in  the  RVCQ3D  code,  identical  test  cases  were 
run  with  the  Baldwin- Lomax  and  RNG  turbulence  models.  Also  a 
fully  laminar  flow  calculation  was  included  for  comparison. 
Unfortunately,  the  RNG  model  did  not  give  a  converged  solution 
after  3000  iterations.  In  Figure  8,  the  convergence  history 
shows  that  the  RNG  turbulence  model  was  just  calculating  the 
laminar  flow  condition. 
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The  flow  velocity  vector  around  the  trailing  edge  for  the 
Baldwin- Lomax  turbulence  model  and  RNG  turbulence  model  are 
shown  in  Figures  9  and  10,  respectively.  The  velocity  vector 
plot  of  the  Baldwin- Lomax  model  solution  shows  a  near- 
symmetric  double  recirculation  region  aft  of  the  trailing  edge 
in  the  turbulent  wake.  The  extent  of  this  recirculation  region 
is  approximately  1.2  trailing  edge  thicknesses  in  the 
streamwise  direction;  however,  the  velocity  vector  plot  of  the 
RNG  model  solution  was  asymmetric.  The  suction- side  boundary 
layer  separated  because  it  was  still  laminar  near  the  trailing 
edge  and  this  resulted  in  the  large  recirculation  region  on 
the  suction  side  which  extended  aft  of  the  trailing  edge  to 
approximately  2.5  trailing  edge  thicknesses  downstream  of  the 
blade.  Separation  started  from  about  1/5  of  the  chord  length 
from  the  trailing  edge  on  the  suction  side  surface.  This  is 
also  the  evident  in  the  laminar  flow  test  case. 
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Figure  9   Velocity  vector  plot  of  B-L  model 
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Figure  10   Velocity  vector  plot  of  RNG  model 
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1.   BLADE  SURFACE  PRESSURE  DISTRIBUTION 

All  of  the  following  comparisons  between  the  B-L  an 
and  RNG  turbulence  model  are  made  for  conditions  after  1000 
iterations,  as  this  is  almost  at  the  best  convergence  of  the 
RNG  model.  However  the  B-L  results  could  be  substantially- 
improved  if  allowed  to  continue  to  full  convergence  beyond 
3000  iterations.  In  Figure  11,  the  prediction  of  the  blade 
surface  static  pressure  to  inlet  total  pressure  ratio  is 
shown. 
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Figure  8  showed  that  the  RNG  turbulence  model  gave  a 
convergence  history  similar  to  the  laminar  flow  solution,  and 
Figure  10  showed  the  laminar  flow  behavior  of  the  RNG  model  at 
the  trailing  edge.  But  Figure  11  shows  that  the  RNG  turbulence 
model  produced  a  different  pressure  distribution  to  the 
laminar  flow  case.  On  examination  of  the  results  over  the 
whole  computational  domain,  the  turbulent  eddy  viscosity 
around  the  airfoil  out  to  the  freestream,  was  found  to  be 
equal  to  zero.  However  the  eddy  viscosity  was  not  equal  to 
zero  in  the  wake  region.  Since  the  mixing  length  formulation 
is  modified  in  the  wake  region,  and  is  different  from  the 
calculation  in  the  region  around  the  airfoil,  this  would 
explain  why  the -RNG  turbulence  solution  was  not  the  same  as 
the  laminar  flow  solution.  And  because  of  this,  turbulent  flow 
was  actually  calculated  in  the  wake  region.  The  RNG  solution 
gave  a  better  prediction  of  the  pressure  distribution, 
particularly  on  the  suction  surface  of  the  blade. 
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In  Figure  13,  the  flow  velocity  to  critical  velocity 
ratio  in  the  circumferential  direction  at  50  percent  chord 
from  the  pressure-side  to  suction-side  is  shown. 
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Figure  13  Turbulence  model  comparison  of  velocity- to- critical 
velocity  ratio  at  50  percent  chord 

The  Baldwin- Lomax  and  RNG  turbulence  model  solutions 
at  the  50  percent  chord  position  are  comparable. 
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2.  MIDCHORD  FLOWFIELD  COMPARISON 

In  Figure  12,  at  the  50  percent  chord  position,  the 
blade-to-blade  distribution  of  Cp  from  pressure-side  to 
suction-side  is  shown.  Here,  all  three  test  cases  show 
comparable  results.  Notice  that  the  RNG  turbulence  model  is 
slightly  different  from  the  other  two  test  cases  in  the  near 
suction- side  region.  This  difference  is  likely  to  be  the 
effect  of  the  recirculation  on  the  suction- side  of  the  blade 
near  to  the  trailing  edge. 
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Figure  12  Turbulence  model  comparison  of  Cp  prediction  at  50 
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For  completeness,  the  predicted  blade-to-blade  flow 
angle  distribution  in  the  circumferential  direction  from  the 
pressure-side  to  suction-side  is  compared  to  the 
experimentally  measured  data,  and  is  shown  in  Figure  14. 
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3.  WAKE  PREDICTION 

The  wake  prediction  with  the  Baldwin- Lomax  and  RNG 
turbulence  models  and  the  laminar  flow  calculation,  are  shown 
in  Figure  15 . 


2 

> 


17.5       20.0  22.5  25.0  27.5  30.0       32.5 

Circumferential  Position,  Theta,  Deg 


Figure  15   Turbulence  model  comparison  of  wake  prediction, 
from  pressure- side  to  suction- side 

The  Baldwin- Lomax  and  RNG  turbulence  models  both  give 
a  good  prediction  of  the  freestream  velocity  magnitude,  but 
the  Baldwin- Lomax  model  gives  the  best  wake  center  line 
position  and  the  velocity  deficit. 
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C.   DISCUSSION 

As  the  RNG  turbulence  model  case  did  not  completely 
converge  after  a  large  number  of  iterations,  two  major 
parameters  in  the  quartic  equation  which  are  inside  the 
Heaviside  function  were  investigated.  These  parameters  were 
the  mixing  length  and  the  shear  stress.  Since  the  RNG  model 
only  solved  laminar  flow  during  the  initial  1000  time  steps, 
this  meant  the  Heaviside  function  inside  the  quartic  equation 
was  always  negative. 

The  constant,  Cc/  in  the  Heaviside  function  was  initially 
set  equal  to  200.  In  order  to  reduce  the  contribution  of  the 
negative  term  in  the  Heaviside- function  calculation,  the 
smallest  value  of  Cc,  75,  was  used. 

Figure  16  shows  the  comparison  of  the  predicted  mixing 
length  for  both  the  Baldwin- Lomax  and  RNG  models.  This  is 
shown  for  the  grid  line  190,  which  is  normal  to  the  blade 
surface,  close  to  the  trailing  edge  (which  is  at  line  200) . 
This  comparison  is  done  at  this  location,  because  it  was  felt 
that  the  boundary  layer  should  be  fully  turbulent  at  this 
station  as  predicted  with  the  Baldwin- Lomax  model.  Transition 
was  determined  to  take  place  at  the  165  grid  line,  which  is 
located  at  approximately  60  percent  chord.  Both  models  reach 
a  maximum  of  Lf  inside  the  boundary  layer  and  these  values  are 
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Figure  16     Comparison   of   B-L   and   RNG   mixing   length 
distributions  normal  to  the  suction  surface 


within  30  percent  of  each  other.  If  we  assume  that  the 
Baldwin- Lomax  model -derived  mixing  length  is  correct  then  the 
RNG  model  mixing  length,  which  is  based  on  the  boundary  layer 
thickness,  6,    seems  to  be  correct. 
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In  Figure  17,  the  shear  stress  distribution  versus  normal 
distance  to  boundary- layer  thickness  ratio  is  shown.  Here  an 
obvious  discrepancy  in  the  level  of  prediction  of  r  is 
evident,  which  most  probably  resulted  in  the  Heaviside 
function  not  going  positive  and  therefore  not  initiating 
turbulent  flow  calculations. 


Figure  17     Comparison   of   B-L   and   RNG   shear   stress 
distributions  normal  to  the  suction  surface 
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An  extra  check  can  be  made  on  the  derivation  of  the  mean 
strain  rate,  equation  (24)  .  The  vorticity  calculated  by  the 
Baldwin- Lomax  turbulence  model  at  the  wall,  and  the  mean 
strain  rate  at  the  wall  around  the  blade  surface,  are  shown 
compared  in  Figure  18 .  They  are  in  agreement  everywhere  except 
at  a  few  points  close  to  both  leading  edge  and  trailing  edge. 
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Figure  18   Shear  strain  and  vorticity  distribution  on  the 
blade  surface 
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The  same  quantities  plotted  along  the  190  grid  line  is 
shown  in  Figure  19.  For  the  near-wall  region,  the  first  15 
grid  points  off  the  blade  surface  are  identical. 
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Figure  19  Shear  strain  rate  and  vorticity  distribution  normal 
to  the  suction  surface 


Since  the  values  of  computed  vorticity  and  mean  strain 
rate  are  exactly  the  same  at  the  blade  surface,  this  verifies 
that  the  derivation  and  coding  of  equation  (24)  are  correct. 
As  the  vorticity  at  the  wall  is  equal  to  the  normal  velocity 
gradient  and  the  mean  strain  rate  is  also  equal  to  this 
quantity. 
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V.   CONCLUSION 

The  annular  turbine  cascade  flowf ield  calculation  with  the 
RVCQ3D  code,  which  is  based  on  the  Baldwin- Lomax  turbulence 
model,  is  grid  dependent.  Increasing  the  number  of  grid  points 
around  the  body  is  a  more  efficient  way  to  improve  the 
solution  than  increasing  the  number  of  points  normal  to  the 
body.  A  quasi- three-dimensional  calculation  which  includes  the 
streamtube  variation  information  is  strongly  recommended, 
particular  for  this  annular  cascade  test  case. 

Initial  indications  were  that  the  RNG  model  showed  similar 
convergence  (a  logical  first  check)  to  the  Baldwin- Lomax 
turbulence  model  for  the  first  1000  iteration,  after  which  the 
RNG  model  diverged.  The  divergence  of  the  RNG  model  could  be 
expected  as  the  model  is  highly  non- linear,  and  most  probably 
less  robust  than  the  Baldwin- Lomax  model.  The  comparison  of 
the  computed  flowfield  at  1000  iteration,  showed  that  the 
Baldwin- Lomax  model  was  overall  better  than  the  RNG  model. 
However/  upon  further  investigation  it  was  determined  that  the 
RNG  model  only  computed  laminar  flow  over  the  blade  surface 
and  turbulent  flow  in  the  wake.  As  the  RNG  model  did  quite 
well  in  the  wake  region,  it  should  be  further  investigated  as 
a  possible  wake  modification  to  the  Baldwin- Lomax  model. 

The  RNG  model  calculation  was  initiated  by  assuming 
laminar  flow  initially  as  described  by  Kirtley  [Ref.  4],  and 
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this  resulted  in  the  Heaviside  function  always  being  negative. 
Another  method,  which  should  be  tried  is  to  begin  from  the 
high- Reynolds  number  limit  as  proposed  by  Lund  [Ref.  3]  and 
use  the  RNG  model  to  update  the  provisional  eddy  viscosity; 
i.e.,  begin  with  a  positive  Heaviside  function. 
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APPENDIX  A    DATA  FOR  GOLDMAN'S  CASCADE 

The  complete  250  X  60  C-type  grid  for  Goldman's  cascade  is 
shown  in  Figure  Al .  Data  for  the_ variation  of  the  streamtube 
thickness  through  Goldman's  cascade  are  given  in  Figure  A2 . 
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Figure  Al  Geometry  of  C-type  250  X  60  grid 
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Figure  A2   The  streamtube  variation  data 
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APPENDIX  B    RESULTS  USING  A  97  X  31  GRID 

Results  are  given  for  calculations  of  the  flow  through 
Goldman's  cascade  using  a  97  X  31  C-type  grid.  The  effect  of 
including  the  effects  of  streamtube  contraction  on  the 
convergence  history  is  shown  in  Figure  Bl,  the  pressure  ratio 
distribution  around  the  blade  is  shown  in  Figure  B2 ,  and  on 
the  calculated  velocity  in  the  blade  wake  is  shown  in  Figure 
B3. 
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Figure  Bl   Effect  of  streamtube  variation  on  convergence 
history  for  the  97  X  31  grid 
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Figure  B2   Effect  of  streamtube  variation  on  static  pressure 
to  inlet  total  pressure  ratio 
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Circumferential  Direction 


Figure  B3  Effect  of  streamtube  variation  on  wake  velocity  to 
critical  velocity  ratio 
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APPENDIX  C    T106  [Ref.  19]  TEST  CASE 

As  another  possible  test  case  for  the  RNG  model,  a  cascade 
of  turbine  rotor  blades  was  considered.  This  section  deals 
with  the  attempt  to  establish  a  baseline  solution  using  the 
Baldwin- Lomax  turbulence  model  of  the  extensive  experimental 
data  set  as  presented  by  Hoheisel  [Ref.  19] . 

The  T106  test  case  is  a  subsonic  turbine  cascade  which 
consisted  of  7  blades  with  a  chord  length  of  c  =  100  mm.  The 
blades  had  an  aspect  ratio  of  3 .  A  set  of  experimental  data 
for  a  Reynolds  number  of  5  x  105  and  0.8  percent  inlet 
turbulence  intensity  was  used  in  this  baseline  calculation. 
The  geometry  of  the  turbine  blade  is  shown  in  Figure  CI,  with 
a  C-type  grid  mesh  of  200  X  30. 

The  RVCQ3D  code  and  an  updated  version  of  RVCQ3D,  with 
consistent  non-dimensionalization  throughout  (by  Tweedt)  ,  were 
both  used  to  compute  the  whole  flow  field.  Baldwin -Lomax 
turbulence  modelling  was  used  in  each  code.  The  boundary 
conditions,  ambient  pressure,  total  "temperature,  pressure 
ratio,  the  leading  edge  flow  velocity,  and  the  inlet  and  exit 
flow  angles  were  provided  in  the  experimental  data  set. 
Unfortunately,  neither  RVCQ3D  nor  Tweedt 's  code  could  satisfy 
the  experimentally  measured  inlet  and  exit  flow  conditions. 
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Figure  CI  Geometry  of  the  subsonic  turbine  cascade,  T106 

The  inlet  velocity  was  varied  over  a  wide  range  of  parameters 
in  the  experiment.  However,  the  RVCQ3D  code  could  not 
reproduce  complete  the  isentropic  velocity  at  the  inlet 
boundary. 
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1.   Cp  COMPARISON 

In  Figure  C2 ,  both  codes  are  seen  to  generate  a  better 
distribution  of  Cp  on  pressure  side  than  on  the  suction  side. 
Overall,  Tweedt's  code  seems  to  give  better  results  than  the 
RVCQ3D  code. 
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Figure  C2   Initial  prediction  of  the  blade  surface  static 
pressure  coefficient  distribution  around  the  blade 
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2.   SKIN  FRICTION  COMPARISON 

The  skin  friction  coefficient  (Cf)  distribution  on  the 
suction  side  of  the  blade  is  shown  in  Figure  C3.  Cf  was 
calculated  from  the  ratio  of  the  shear  stress  to  the  local 
dynamic  pressure  at  the  wall.  Neither  prediction  approaches 
the  experimental  data. 
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